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Introduction 
The introduction to our Fall 2010 ELEC 301 project on surround sound 
simulation. 


Introduction 


Background 


In today’s rapidly growing entertainment industry, many people desire a 
greater degree of emersion whether it be through 3-D visuals popping out of 
the screen or surround sound audio making them feel as though they are at 
the heart of the excitement. Unfortunately, the means to implement these 
systems can often be out of reach for individuals due to the cost of such 
systems or their size. Surround sound systems, while not as expensive as 
one of the new 3-D TVs, are still a significant investment both in money 
and the time needed to correctly set up the system. Additionally, these 
systems are by no means portable. 


Project Description 


Through our project we sought to address both the cost and size issue of 
surround sound systems by using in-ear microphones to sample a test signal 
from which a person’s particular head-related transfer functions (HRTFs) 
could be calculated. These transfer fnctions could then be applied to input 
signals and the resulting output signals, when played through stereo 
headphones, would give the appearance of sounds coming from various 
directions in a circle around the person’s head. 


Test Sampling and Channel Characterization 
The means by which we characterized our speaker to mic channel using a 
test signal. 


Sample Signal Generation 


For our sample signal we decided on a chirp covering the range of human 
hearing, 15 Hz to 20000 Hz. This signal was created using the Matlab chirp 
command as follows: 


%% Chirp Function 


Fs = 44100; 

fO = 15; 

fi = 20000; 

t = 0:1/Fs:1-1/Fs; 

y = chirp(t,f0,1, 1); 


where f0 and f1 are the edges of the spectrum of human hearing. The vector 
t has 44100 entries, each separated by 1/44100. Taken as a time vector and 
with the end time of the chirp set as 1 second, this creates a chirp which 
goes from 15 Hz to 20000 Hz in 1 second at intervals of 1/44100 seconds, 
which allows for a clean transition from 15 Hz to 20000 Hz at the standard 
sampling rate of 44100 Hz. 


Finding the Transfer Function of the Speaker to Mic Channel 


Microphone and Speaker Equipment used for signal sampling 


In order to find HRTFs we first had to know the transfer function of the 
speaker to mic channel through air. We did this by first playing our chirp 
signal at a mic held approximately 3 feet away and recording the mic 
output. We did this twice and averaged the two outputs to remove any 
anomalies. Then using the property that for an LTT filter: 

Equation: 


channel (t) * in(t) = out(t) 
or in terms of frequency domain: 
Equation: 


H channel (w) . IN(w) = OUT(w) 


Equation: 


A Gis (w) = OUT(w) /IN (w) 


where AH channel(w)is the transfer function of the system and IN(w) and 
OUT(w) are the fourier transforms of the input and output signals 
respectively. Thus, by taking the fft of the mic output and pointwise 
dividing by the fft of the chirp signal (3), we found the transfer function of 
the speaker to mic channel. 
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Finding Directional HRTFs 
Finding directional transfer functions by means of in-ear signal sampling. 


Calculating HRTFs for Signals from Various Directions 


Directional Signal Sampling 


In order to find directional HRTFs for our test subject, we first had to make 
recordings of the signals as heard by the test subject coming from different 
directions. In order to do this, we placed our microphone in the ear to be 
tested, as shown below. 


Custom, In-Ear 
Microphone 


In-ear microphone setup for the left ear 


We then played our chirp signal at the test subject from various directions 
with the speaker at the same distance as in the initial channel 


characterization. As in the initial channel characterization, we sampled 
twice for each direction and averaged the two results. After sampling for 
that particular ear, we then switched the microphone to the other ear and 
sampled the same directions for the other ear. 
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The five sampling directions used to find HRTFs 


HRTF Calculation 


In the initial channel characterization, the transfer function was 
the only thing acting on the chirp signal. Now, with the microphone situated 
in the test subject’s ear, both and the directional HRTF for that 
particular ear act on the chirp signal. In other words: 


Equation: 


Since we already calculated and during the initial 
channel characterization, and we can find for each particular 
ear/direction combination by taking the fft of our recorded outputs, we can 
therefore calculate using the equation: 

Equation: 


Surround Sound Simulation 
Generating the illusion of directional sound using calculated directional 
transfer functions. 


Directional Test Signal Creation 


In order to test the accuracy of our calculated HRTFs, we used 5.1 surround 
sound test files used to calibrate home theater systems. We first had to 
downsample these files from their original sampling rate of 48000 Hz to our 
sampling rate of 41000 Hz. Since the files were all meant to be played at 
the same time, each file only played sound part of the time. Therefore, we 
had to cut out the unused parts of each file and create new sound files which 
included only the testing segment of each directional file. After this process 
was complete, we were left with five signals, each of which specified which 
direction the sound was supposed to come from. 


Directional Sound Playback 


We then simulated these signals coming from particular directions using the 
equation from the HRTF calculation: 


Equation: 

A 8) w HH WwW W 
In this case H is the earbud to ear channel transfer function. 
Unfortunately, we weren’t able to characterize this channel with our 
available resources, so in this instance we assumed that H . This 
left us with the following equation: 
Equation: 

w H W 7) 

By applying this equation to each calculated H w for both ears we 


were able to generate a series of synthesized sound files for each ear which 


when played in stereo would mimic a cound coming from a particular 
direction by using the following code: 


EZOUt= (Sy Nie SYNPIBECIE, -SyYNinZeReE, 
Sy Nitin RR es oY Niet ee 


L_OUT=L_OUT'; 


ROOUT=(SYNTHUFESR, SYNTHOFC R, SYNUH_ERER, 
SYNTnoRRIRe voy Nino RE R | 


R_OUT=R_OUT'; 
FINAL_OUT=[L_OUT, R_OUT]; 
sound(FINAL_OUT, Fs); 


In this code, the SYNTH terms are the vectors corresponding to the filtered 
test signals for front-left, front-center, front-right, rear-right, and rear-left 
directions. With the left and right ear outputs arranged as shown in 
FINAL_OUT, the sound command plays the first set of terms to the left 
speaker/headphone and the second set of terms to the right 
speaker/headphone as we desired. 


Results, Conclusions, and Future Applications 
The results, conclusions, and future applications of this surround sound 
simulation project. 


Project Wrap-up 


Results 


When the synthesized sound playing command is executed, the results are 
encouraging. Obviously, the surround sound simulation works best for the 
person who was used in the HRTF calculation step as each person will have 
different HRTFs. The illusion of sounds coming from different directions is 
clearly discernible to other listeners, but generally the directions are neither 
as distinct nor as accurate. The resulting sound is also somewhat noisy, due 
to the testing environment and the quality of the equipment used. 


Conclusion 


While this method of surround sound simulation wasn’t quite as convincing 
as a home theater setup might have been, it was quite accurate considering 
the number of HRTFs used, the testing environment, and its’ ability to work 
in part for more people than simply the HRTF test subject. 


Future Directions/Applications 


This simulation method, if perfected, could be readily scaled to calculate 
more directional transfer functions. This would allow the simulation of 
sound with greater directional fidelity. Possible means of improving this 
method include: 


e Recording in a less noisy, anechoic chamber 

e Using higher quality sound production and recording devices with 
more linear transfer functions 

e Recording at higher sampling rates 


e Measuring the transfer function between each headphone ear bud and 
the user’s inner ear, and incorporating this into the directional transfer 
function 

e More accurately measuring speaker location/angle to better simulate a 
home theatre environment 

e Use of a nonlinear chirp signal to account for the nonlinearity of 
human hearing 


An interesting hypothetical application of this technology involves noise- 
cancelling headphones. Such headphones have a small microphone 
embedded in them for the purpose of sampling environmental noise. These 
microphones could conceivably be used to sample a particular user’s 
HRTFs if an external speaker was available. These, along with software, 
could then be used to simulate a multi-directional speaker configuration. 
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